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We develop a theory of the effective disorder temperature in glass-forming materials driven away 
from thermodynamic equilibrium by external forces. Our basic premise is that the slow config- 
urational degrees of freedom of such materials are weakly coupled to the fast kinetic/vibrational 
degrees of freedom, and therefore that these two subsystems can be described by different tem- 
peratures during deformation. We use results from the preceding paper on the nonequilibrium 
thermodynamics of systems with internal degrees of freedom to derive an equation of motion for the 
effective temperature and to learn how this temperature couples to the dynamics of the system as 
a whole. 



I. INTRODUCTION 

The effective disorder temperature is emerging as an 
essential ingredient in theories of nonequilibrium phe- 
nomena in amorphous materials [j], 0, y|, 0, [![• Much 
of our own most recent work is based on papers by Liu 
and colleagues, especially [a0]- The preceding authors 
showed that, in systems driven by external forces, slowly 
relaxing degrees of freedom may fall out of equilibrium 
with the thermal reservoir and, as a result, the temper- 
ature associated with those slow degrees of freedom dif- 
fers from the reservoir temperature. More recently, the 
effective temperature has been used extensively in shear- 
transformation-zone (STZ) theories of large-scale plastic 
deformation in molecular glasses [j| [tj [T(J, [111, [l2|, EH- 
Perhaps its most remarkable success along these lines has 
been in explaining the nature of shear-banding instabili- 
ties in such systems [14], [l5| . 

Despite these successes, it has never been clear - at 
least not to the present authors - whether the effective 
temperature is truly a "temperature" in the thermody- 
namic sense, or whether it merely resembles a temper- 
ature in some aspects but differs in others. We have 
argued, e.g. in [{j, M, H3, [l3| . that a noncrystalline mate- 
rial should be characterized in large part by its internal 
state of disorder, which we identify with the configura- 
tional entropy. Then the intensive quantity thermody- 
namically conjugate to this entropy, i.e. the derivative of 
the configurational energy with respect to the configura- 
tional entropy, must have the properties of a temperature 
in the sense that systems with higher temperatures are 
more highly disordered, and neighboring subsystems are 
in equilibrium with each other if they are both at the 
same temperature. 

Our goal here is to make these ideas more precise. By 
doing so, we wish to determine what roles the effective 
temperature can play in dynamical theories of glassy phe- 
nomena, and what its limitations might be. Questions to 
be addressed include: Under what circumstances is the 
effective temperature different from the ordinary reser- 



voir temperature? What are the forms of the first and 
second laws of thermodynamics in systems where both or- 
dinary and effective temperatures are well defined? What 
is the equation of motion for the effective temperature? 
For the most part, our answers confirm the guesses made 
in earlier publications, but there are some differences. 

To answer these questions, we have had to reexamine 
some basic principles of thermodynamics, especially the 
meaning of the second law in a non-standard situation 
of the kind we encounter here; and we also have had to 
understand the thermodynamic roles played by internal 
degrees of freedom in nonequilibrium situations. We have 
reported progress along those lines in the preceding paper 
|16| . Much of the present analysis is based on the results 
described there. 

In Sec. [TT] of this paper, we introduce the hypothesis 
that a glassy system consists of two weakly coupled sub- 
systems, one of which contains the slow configurational 
degrees of freedom that are characterized by the effective 
temperature. In order to focus as sharply as possible on 
the basic thermomechanical issues, as opposed to ques- 
tions about spatial variations or time dependent orienta- 
tions of the stress, we specialize to the case of a spatially 
uniform system. We make several other simplifying as- 
sumptions, to be specified later, for the same reason. We 
derive the the two-temperature equations of motion in 
Sec. IIII1 and conclude in Sec. IIVI with remarks about 
the theory. 



II. SEPARABLE SUBSYSTEMS 

Our basic premise is that the degrees of freedom of a 
classical (h = 0) noncrystalline material, either below or 
not too far above a glass transition, can to a good ap- 
proximation be separated into two parts: first, a set of 
configurational coordinates describing the mechanically 
stable positions of the molecules in their inherent states 
[l7l [l8l [T^| . and, second, the remaining kinetic and vibra- 
tional variables. We assume that the two sets of degrees 
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of freedom describe weakly coupled subsystems of the 
system as a whole, in analogy to the usual thermody- 
namic analysis in which two neighboring subsystems of a 
larger system are brought into weak contact. In the latter 
situation, if the temperatures of the neighboring subsys- 
tems differ from each other, we know from the second law 
of thermodynamics that heat flows from the hotter to the 
colder subsystem; and we can compute the rate of heat 
flow either phenomenologically by defining a transport 
coefficient, or sometimes by computing that coefficient 
from first principles. We propose to repeat such an anal- 
ysis for configurational and kinetic/vibrational degrees of 
freedom in weak contact with each other. 

The configurational and kinetic/vibrational subsys- 
tems are distinguished by their different time scales. The 
kinetic/vibrational degrees of freedom are fast; they move 
at molecular speeds. The configurational degrees of free- 
dom are slow; they change on time scales no shorter than 
those associated with molecular rearrangements in glassy 
materials. If these two time scales are comparable with 
one another, then the system is liquidlike, and the con- 
cept of an effective disorder temperature is irrelevant; 
but, as will be seen, the separation between fast and slow 
variables may remain useful. 

To be more specific about the separability of these sub- 
systems, consider a potential-energy landscape in Nd- 
dimensional configuration space, where N is the number 
of molecules and d is the dimensionality. The inherent 
states are the molecular configurations at the local en- 
ergy minima in this space; and the configurational sub- 
system, apart from uniform translations and rotations, is 
defined by the denumerable set of these states. A molec- 
ular rearrangement of the kind that occurs during plastic 
deformation corresponds to the motion of a system point 
from the near neighborhood of one inherent state to that 
of another. 

If the total energy of the system lies below all but a 
small fraction of the energy minima, then the intersec- 
tion of the potential-energy landscape with the surface of 
constant energy consists primarily of isolated, compact, 
manifolds enclosing the low- lying inherent states. Very 
few energy-conserving, deterministic trajectories connect 
the neighborhoods of different inherent states. For a 
large enough system that is not totally jammed, however, 
there always must exist some such trajectories; neigh- 
boring molecules always can move around each other if 
a large enough number of the other molecules move out 
of their way, as happens in the excitation-chain mecha- 
nism for example [201 ] . These connecting trajectories - 
narrow filaments in configuration space - occupy only 
a small fraction of the constant-energy manifold. Thus, 
in physical space, they correspond to highly improbable, 
collective events. The probabilities of such events are 
enhanced when the transitions are assisted by thermal 
fluctuations, which are not yet part of this picture but 
will enter shortly. 

By definition, the motions of the kinetic/vibrational 
subsystem occur only on trajectories within the compact 



manifolds and not on the filamentary connections be- 
tween them. Our principal assumption is that the latter 
motions can accurately be approximated by weak, ther- 
modynamically constrained, couplings between the con- 
figurational and kinetic/vibrational subsystems. That is, 
molecular rearrangements in the configurational subsys- 
tem - transitions between inherent states - are driven 
by noise generated in the kinetic/vibrational subsystem; 
conversely, heat generated during irreversible rearrange- 
ments flows through the kinetic/vibrational subsystem to 
a thermal reservoir. 

We further assume that we can accurately approximate 
the total potential energy near any one of the minima by 
the inherent-state energy plus the harmonic energy of 
small-amplitude molecular vibrations about that state. 
Then, for a system of N molecules in the neighborhood 
of the i/'th inherent state, the total energy H is given by 

H = H c {r„} + H K {p, 5r v }, (2.1) 

where 

H K {p,8r v }=K{p} + W{5r v }. (2.2) 

Here, K {p} is the kinetic energy as a function of the set 
of N momenta {p}, -ffc{ T v} is the inherent-state energy 
at the i/'th minimum located at {r v }, and ^{Jr,,} is 
a quadratic form in the molecular displacements {Sr u } 
measured from that minimum }21| . Note that the sepa- 
ration of the subsystems postulated in Eq. (|2.ip is math- 
ematically exact if the vibrational spectrum associated 
with Hk {p, 8r v } were independent of the inherent-state 
index v\ this is not quite true in general, but is a rather 
good approximation below the glass transition in the ab- 
sence of irreversible volume deformation • This issue 
will be discussed in detail in Sec. Mil below. 

The configurational part of the internal energy, He, 
must depend on the size and shape of the system; i.e. He 
must contain at least the macroscopic part of the elastic 
energy. Since, by construction, both subsystems occupy 
the same region of physical space, we must exclude the 
macroscopic displacement modes from Hk- Doing so is 
consistent with the separation of time scales because elas- 
tic modes that span the whole system have frequencies 
that are vanishingly small in the limit N — > oo. We 
thus assume that the bulk elasticity resides primarily in 
the configurational subsystem. The kinetic/vibrational 
subsystem consists only of the high-frequency, localized 
modes that, in a disordered noncrystalline material, do 
not probe the boundaries of the system and are sensi- 
tive only to the average spacing between the molecules. 
Therefore, we assume that the kinetic/vibrational sub- 
system has an elastic bulk modulus and contributes a 
partial pressure to the system as a whole; but it has no 
shear modulus. 

Note that we are maintaining the analogy between the 
present picture of weakly coupled configurational and ki- 
netic/vibrational subsystems, and the usual thermody- 
namic picture of spatially neighboring subsystems. In the 
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usual picture, the neighboring dynamical systems have 
coordinates in common only for the molecules that inter- 
act with each other across the surface of contact, and the 
coupling is weak because it involves only a small fraction 
of the total number of degrees of freedom. Here, because 
we assume that the basins of attraction of the inherent 
states do not appreciably overlap or, equivalently, that 
the vibrational amplitudes are small, the coupling is weak 
because only a small fraction of the configuration space 
contains trajectories that link different inherent-states. 



III. TWO-TEMPERATURE 
THERMODYNAMICS 

A. Internal Energy and Entropy 

In accord with the preceding discussion, we consider 
the deformation of a spatially uniform, amorphous sys- 
tem, in contact with a thermal reservoir at tempera- 
ture Or. We assume that Or is either below or not 
too far above the glass temperature g , so that the two- 
temperature decoupling can take place. As in we 
express temperatures in units of energy, and set Boltz- 
mann's constant ks equal to unity. 

Our goal of identifying temperature-like quantities re- 
quires that we start the statistical analysis by working in 
a microcanonical ensemble where the internal energy of 
each subsystem is a function of its entropy - in princi- 
ple, an independently computable quantity obtained by 
counting states. Thus we write the total, extensive, inter- 
nal energy of the system, including a thermal reservoir, 
in the form 

U tot U c (S c ,E el , {A a }) + U K {.S K , V*) + U R (S R ) , 

(3.1) 

where Uc and Uk, respectively, are the configurational 
and kinetic/ vibrational internal energies obtained from 
the energy functions He and Hk discussed in Sec. HH Sc 
and Sk are the respective entropies, and {A Q } denotes 
a set of internal state variables. Variations of the A a 
describe irreversible changes in the state of the system 
which may be - but are not necessarily - coupled directly 
to shape deformations. Ur is the energy of the thermal 
reservoir, which we assume to be strongly coupled to the 
kinetic/ vibrational subsystem. 

E e i is a symmetric elastic strain tensor that, in princi- 
ple, includes both deviatoric (shear) and volumetric de- 
formations. We maintain that only the reversible part of 
the strain can appear as an independent argument of the 
internal energy. The inelastic deformation can play no 
role here, because it must be defined as the displacement 
from some initial reference configuration and therefore is 
determined by the entire history of prior deformations. 
We insist that the motion of the system at any given 
time be determined entirely the external forces acting 
on it and the current values of the internal state vari- 
ables. In other words, the memory of prior deformations 



is carried entirely by those state variables. The inelastic 
deformation itself can be determined only by integrating 
the equations of motion for those quantities. 

In using the elastic part of the deformation as an inde- 
pendent argument of the internal energy, we are assuming 
that an elastic strain tensor can be identified unambigu- 
ously This is a controversial assumption, for example, 
see the discussion in [23| ; but it is not crucial for the 
main theme of this paper. Therefore, for simplicity, we 
restrict ourselves to situations in which elastic displace- 
ments are small and linear elasticity is accurate; but we 
consider arbitrarily large inelastic deformations. 

As argued in Sec. [TT1 the kinetic/vibrational subsystem 
has a bulk modulus, and therefore the internal energy 
Uk in Eq. (|3.ip can depend on the elastic part of the 
volume deformation, i.e. V e i — V trE e f, but it does not 
depend on E e i itself. The absence of a shear modulus for 
this subsystem follows rigorously from our assumption of 
linear elasticity. To see this, invert Uk(Sk, Vei) to obtain 
the kinetic/vibrational entropy Sk in the form 



S K = S K (U K ,V e i) 



(3.2) 



In the harmonic approximation for Hk in Eq. (|2.1[) . we 

have 



~>K 



3 N 



1+ln 



Uk 
3N 



3N 



^ln(w„) 



(3.3) 



where the u> n are the frequencies of the vibrational 
modes, and we have assumed there to be 3A^ kinetic 
and 3A*" vibrational degrees of freedom. The dependence 
of Sk on the deformation resides in the frequencies u) n , 
which are the eigenvalues of the dynamical matrix - the 
second derivatives of V7{(5tv} in Eq. (|2.2[) . For linear 
elasticity, the u) n are linear in that strain, and there- 
fore the first-order shifts in the normal-mode frequencies 
must change sign when the shear strain changes sign. On 
the other hand, rotational and inversion symmetries in a 
statistically isotropic system preclude such sign changes. 
Thus, the Lu n must be independent of the shear strain to 
first order, and the shear modulus must vanish. 

The preceding discussion implies that the normal- 
mode frequencies oj n depend on the total volume V, and 
not just on its elastic part. Thus, in principle, Uk de- 
pends on the internal state variables. This mechanical 
coupling between the kinetic/vibrational and configu- 
rational subsystems could be accounted for within the 
present formulation, but only at the expense of unneces- 
sary complication. We therefore assume that the depen- 
dence on the set of internal state variables {A Q } resides 
exclusively in the configurational internal energy Uc- 

In analogy to the expression for the internal energy in 
Eq. (|3.ip . the total entropy is given by 

S tot ~ S C (U C , E e i , {A Q }) + S K (U K , V e i) + S R {U R ) . 

(3.4) 

Note that the "nonequilibrium" subscript used in [l6| is 
no longer needed here. We assume that we can invert the 
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expression for any one of these three entropies to obtain 
the corresponding internal energy function in Eq. (|3.1|) 

or vice versa. 

In principle, we have three separate temperatures: the 
effective temperature \i the kinetic/ vibrational temper- 
ature 0, and the reservoir temperature Or = 0. These 
are given respectively by 



X 



fdU c 



\dS c 



'dU 



K 



E e l,{A a } 



\dS K 



OUr 
dS R 



(3.5) 



Note that we depart here from the notation used in [13| 
and other earlier STZ papers, where x was defined to be 
a dimcnsionless ratio of the effective temperature to an 
STZ formation energy in units of Boltzmann's constant. 
Here, % has the units of energy. The thermodynamic 
stresses are 



VT C = 



dUr 



dE e 



1/ ScfAoJ 



PK 



fdU 



K 



\BVel 



(3.6) 



B. The First Law of Thermodynamics 

The first law of thermodynamics for this system is 

VT:D = U tot , (3.7) 

where the total stress T is the sum of the partial stresses, 
T = Tc + Tk- D is the total rate of deformation 
tensor, including both elastic and inelastic parts. Note 
that all the stress quantities can be expressed as a sum 
T = —p 1 + s of hydrostatic (p) and deviatoric (s) con- 
tributions. Breaking Eq. (|3.7|) into its component parts, 
we obtain 

(t/ c - VT C : D) + (U K -VT K : D) + U R = 0. (3.8) 

To evaluate these terms, look first at the configura- 
tional subsystem, and recall that changing the shape of 
the system at fixed values of Sc and {A Q } is by definition 
an "elastic" deformation; there is no irreversible change 
in the internal state of the system. Using this definition, 
we have 



dU c 
dt 



VT C : D 



el * 



(3.9) 



S C ,{A Q } 



where D e i is the elastic part of the rate of deformation 
tensor. Because we assume that the total rate of defor- 
mation is a sum of elastic and inelastic parts, 



D = D e i + D in , 
the configurational component of Eq. 



(3.10) 



VT c :D = X Sc-VT c :D, 



E 



dU c 



(3.11) 



Next, consider the kinetic/ vibrational component of 
Eq. ([H]), where 



fdU 



K 



V dt 



-PkVci 



(3.12) 



Therefore, 
U K - VT k :D 



eS K -Vs K :D dev + PK V m , (3.13) 



is 

the total rate 

V = Vet + V ln 

Sx in the last 
viatoric stress, 
subsystem can 
appear in Eq. 
flow, i.e. when 
Recombinin^ 
(13. 7[) , becomes 



the deviatoric (non-hydrostatic) part of 
of deformation D, and we have used 

to eliminate V e i. Note the appearance of 
equation. If it has no thermodynamic de- 

cf. Eq. (|3.6p . then the kinetic /vibrational 

support a viscous shear stress, which can 
(|3 . 1 3[) in the presence of nonzero shear 

D dev ^Q. 

terms, we find that the first law, Eq. 



xSc 



Using 



1 S K + U R - V T c : An + PKVir 



(3.14) 



V s K : D dev + 



fdU c 



\dA a 



Aa =0 



S a ,E e 



-VT C : D m = -Vs c : D?™+ Pc V m , (3.15) 
we rewrite Eq. (|3.14[) in the form 

X S c + OSk + Ur-V s c : Df™ + pV in (3.16) 



V s K : D 



dev 



E 



dUc 
0A o 



A r 







So, Eel 



Note that the partial deviatoric stress sc couples to the 
deviatoric, inelastic rate of deformation Df^ v , the total 
pressure p couples to the inelastic part of the volume rate 
of change U; n , and the partial deviatoric stress sk couples 
to the total deviatoric rate of deformation D dev . These 
features will reappear below in dissipation inequalities 
based on the second law of thermodynamics. 



C. The Second Law of Thermodynamics 

The second law of thermodynamics requires that the 
total entropy be a non-decreasing function of time 



where 



S c + S K + S R > 0, 



Ur — Or Sr. 



(3.17) 



(3.18) 



Use Eq. (|3.16[) to eliminate Sc in Eq. (|3.17[) . and rear- 
range the terms as follows 
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Vs K :£>* w +W (s 0l p,{A a }) >0, (3.19) 



where 



E 



9A C 



(3.20) 



where W is the analog of the quantity defined in Eq. 
(4.11) of [Hj]. Note that W c is the difference between 
the rate at which inelastic work is being done on the 
configurational subsystem and the rate at which energy 
is stored in the internal degrees of freedom. As in jl6j], it 
is an energy dissipation rate. 

The left-hand side of the inequality in Eq. (|3.19p is the 
sum of independently variable quantities that must sep- 
arately be non-negative if the inequality is to be satisfied 
for all possible motions of the system. Thus, in the spirit 
of Coleman and Noll [HI , we require that 



1 



U R >0 



and 



, s K + — u R ) >o 

Or 



s K : D dev > 



Wc(sc,P,{A Q }) >0 



(3.21) 



(3.22) 



(3.23) 



(3.24) 



The first two of these inequalities pertain to the heat 
flow between the subsystems. We have arranged them 
in a way that reflects the fact that the thermal coupling 
between the configurational and kinetic/ vibrational sub- 
systems is expected to be weak, and the coupling be- 
tween the kinetic/ vibrational subsystem and the reser- 
voir should be strong. To satisfy them, we write 



>S K + — U R 
Or 



Mx,e) i-£ 



x 



3(6,6! 



1 



6r 



(3.25) 



(3.26) 



where A and B are non-negative thermal transport co- 
efficients. Using the definition of the kinetic/ vibrational 
(extensive) heat capacity at constant volume 



>S K = C$t 



(3.27) 



and eliminating Ur between Eqs. (|3.25[) and (|3.26[) . 
we obtain the following heat equation for the ki- 
netic/vibrational subsystem 



C$ 6 = B( 



i) 



l -T R 



x 



(3.28) 



The first term on the left hand side describes heat trans- 
fer between the kinetic/ vibrational subsystem and the 
thermal reservoir, and the second term describes heat 
transfer between the kinetic/ vibrational subsystem and 
configurational subsystem. Our assumptions imply that 
A(x,6) is small because it contains the physics of the 
weak coupling described in Sec. HU On the other hand, 
a large value of B implies that 6 — > 6r such that 



Q R = lim B(6,6 R ) 

B— >oo 



l 



(3.29) 



is finite. In this way, the thermal reservoir controls the 
temperature of the kinetic/vibrational subsystem. 

The third term on the left hand side of Eq. (f3~19|) 
must, by itself, be non-negative for arbitrary values of the 
deviatoric part of the rate of deformation tensor, D dev . 
We satisfy this requirement by assuming a linear relation 
of the form 



s K = r\D 



dev 



(3.30) 



where the shear viscosity r\ is a scalar for our isotropic 
material. 

The fourth inequality, Eq. (13.24j) . pertains to the 
irreversible dynamics of the configurational subsystem. 
This inequality, which is a generalization of the standard 
Clausius-Duhem inequality 25] to the two-subsystems 
situation, must be satisfied by any constitutive model. 
The implications of this inequality for the STZ theory of 
amorphous plasticity [IJl, [TH [U, [H, HH are discussed 
in the following paper [261 ]. 



D. Equation of Motion for x 

The crux of our two-temperature analysis is an equa- 
tion of motion for the effective temperature x- This is 
a heat-balance equation, which we obtain by rewriting 
the first law of thermodynamics as stated in Eq. (|3.16D . 
Define the effective (extensive) heat capacity at constant 
volume 



X$C 



W X 



(3.31) 



Then use this definition and Eq. ([3~2"S"]) in Eq. $5J§§ to 
obtain 

Cfx = Wc(sc,P,{^ a }) 

+V V D dev : D dev +A{ Xl 6) (l - |) . (3.32) 

Here 

-A( X) 6) (l-|) ^6S k + Ur = -Q (3.33) 

is the rate at which heat is flowing to the ki- 
netic/vibrational subsystem and the thermal reservoir. 
The approximation in Eq. (|3.33p is accurate when Ok Or 
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and the kinetic/ vibrational subsystem serves as the ther- 
mal reservoir. Moreover, if the thermal reservoir has a 
large heat capacity (as is commonly assumed), Eq. (|3.28|) 
implies 

9 = 0, Q R = Q . (3.34) 

Under these circumstances, Eq. (|3.32| is a complete de- 
scription of the energy flow in this system. 

In the absence of deformation, i.e. when the term pro- 
portional to A(x, 0) is all that remains on the right-hand 
side of Eq. (|3.32p , or when A(x,9) is large, then x — ► 9 
and our two-temperature formulation reverts to an or- 
dinary one-temperature theory. A single heat equation 
for x = can be obtained by summing Eqs. (|3.28[) and 
(|3.32|) . and identifying Cy = Cy + Cy f as the total heat 
capacity. For example, using this procedure, we recover 
the one-temperature equations of motion for viscoelastic 
volume deformations derived in [la] . 

We reiterate the physical meaning of the different 
terms in Eq. (|3 . 32|) . C e y x 1S the rate of change of the 
heat of disorder. The first term on the right-hand side, 
Wc(s<7)P> {A Q }), is the non-negative rate of energy dis- 
sipation in the configurational subsystem, which we find 
to be the difference between the rate at which inelastic 
work is being done, (Vsc ■ Df^ v —pVi n ), and the rate at 
which energy is being stored in the configurational de- 
grees of freedom {A Q }. The term V T) D dev : D dev is the 
rate of viscous heating in the kinetic/ vibrational subsys- 
tem. Finally, — A(\, 9)(1 — x/&) is the rate at which heat 
flows out of the configurational subsystem. 



can be described by internal state variables. We stress 
that we have provided only a general thermodynamic 
description, and not yet a detailed physical theory of 
such systems. Although we have made several inessen- 
tial simplifying assumptions - spatial homogeneity and 
a simple distinction between clastic and inelastic defor- 
mations, for example - our only fundamental physical 
assumption is weak coupling between the configurational 
and kinetic/vibrational subsystems. 

To make further progress, we must make more spe- 
cific physical assumptions. In particular, we must 
specify physical ingredients of the dissipation function 
Wc(sc,P, {A Q }) defined in Eq. ([3~2U|l . and the thermal 
coupling coefficient A(x,9) defined in Eq. (|3.25p ; and 
we must make use of the second-law constraints in both 
cases. As was shown for a special situation in [la ], non- 
negativity of Wc(sc,P, {A a }) is a generalization of the 
Clausius-Duhem inequality. Inequalities of this kind are 
usually satisfied [25| by assuming the existence of vari- 
ational principles such that the internal variables {A Q } 
move downhill in some free energy landscape. We shall 
see in [26j that, when we interpret the internal state vari- 
ables {A a } to be the STZ dynamical variables, we no 
longer recover a simple variational form of this inequal- 
ity, and that the dynamical structure of the theory is 
much richer than it was in [l6[. Similarly, much of the 
physics of the STZ theory resides in the coupling coeffi- 
cient A{xi 9). It is here where the assumption of weak 
coupling between the subsystems must be described in 
more detail, albeit still in a phenomenological manner. 



IV. CONCLUDING REMARKS 

The main goal of this paper has been to derive the con- 
straints imposed by the first and second laws of thermo- 
dynamics on driven, glass-forming systems whose states 
of configurational disorder can be characterized by an ef- 
fective temperature, and whose nonequilibrium dynamics 
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